Homological characterizations of spiral defect chaos in Rayleigh-Benard convection 
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We use a quantitative topological characterization of complex dynamics to measure geometric 
structures. This approach is used to analyze the weakly turbulent state of spiral defect chaos in 
experiments on Rayleigh-Benard convection. Different attractors of spiral defect chaos are distin- 
guished by their homology. The technique reveals pattern asymmetries that are not revealed using 
statistical measures. In addition we observe global stochastic ergodicity for system parameter values 
where locally chaotic dynamics has been observed previously. 



Characterization of geometric structures (patterns) in 
physical systems can give insight into underlying dynam- 
ics. For simple cases, pattern characterization is easily 
done with a few numbers (e.g. the lattice constants de- 
scribing problems with crystalline symmetries). However 
encoding patterns becomes more difficult as patterns be- 
come more complex. Statistical approaches can be used 
to describe complex patterns in dynamically evolving sys- 
tems [Tl[2]; however, a general methodology for extract- 
ing geometric signatures from complex patterns has been 
lacking. 

We describe an approach that uses algebraic topol- 
ogy to extract consistent and robust geometric proper- 
ties from experimental observations of spatiotemporal 
complexity. We employ homology theory since it is the 
most computable algebraic topological invariant (i.e., re- 
mains constant under deformations and is independent of 
any particular metric) that still provides detailed geomet- 
ric information. Hydrodynamic systems readily produce 
complex spatiotemporal patterns [3 , which are ideal for 
testing these techniques. We use homology to obtain 
a limited set of numbers (called Betti numbers defined 
below) that characterize dynamically relevant features 
in a weakly turbulent state of Rayleigh-Benard convec- 
tion (RBC) known as spiral defect chaos (SDC) [4 (see 
Fig. [T]). In particular, we find that: (a) the geometric 
structures of asymptotic states of SDC are readily dis- 
tinguishable as a function of a control parameter; (b) the 
entropy of the SDC attractors quantifies the evolution 
of the system through different geometric configurations; 
(c) the asymptotic global geometric configurations evolve 
stochastically as opposed to chaotically as identified pre- 
viously [5 , 6 ; (d) symmetry breaking in the flow patterns 
is readily detectable in cases where simple statistical mea- 
sures (e.g., measurements of mean pattern volumes) are 
insensitive to pattern asymmetries; and (e) the boundary 
and bulk patterns can be distinguished. 

We measure convective flow in a 0.069 cm deep hori- 
zontal layer of CO2 gas at a pressure of 3.2 MPa bounded 
above by a 5 cm thick sapphire window and below by 
1 cm thick gold plated aluminum mirror. The lateral 
walls are circular, formed out of an annular stack of filter 




FIG. 1: Shadowgraph images of spiral defect chaos convection 
patterns at reduced Rayleigh numbers e ^ 1.0 (a) and e ^ 
2.0 (b) are converted to binary valued images [(c) and (d), 
respectively] by thresholding the data at the median value of 
intensity. In (c) and (d), the black regions correspond to hot 
upflow, and the white regions to cold downflow. 



paper sheets 3.8 cm in diameter. An electrical resistive 
heater is used to heat the bottom mirror, and the top 
window is cooled by circulating chilled water at 21.2 ^C. 
When the vertical temperature difference, T, across the 
gas exceeds a critical temperature difference, Tc ~ 4 ^C, 
the onset of fluid motion occurs. The flow organizes into 
a pattern of convective rolls where hot regions of the 
gas move upward and cold regions flow downward. The 
convection rolls become spatially disordered and exhibit 
complex time dependence when the system control pa- 
rameter e = (T — Tc)/Tc (the reduced Rayleigh number) 
is sufficiently large. For the present study, the patterns 
are observed for the order of lO^r^^ (r^^ ~ 2.1 seconds 



is the vertical thermal diffusion time) at e ~ 1.0 and 

Shadowgraph visualization of the convecting flow 
yields intensity images, which are used to compute the 
topology of the rolls in SDC (Fig. [T]). The shadowgraph 
technique is sensitive to the refractive index variations in 
the gas, which represents the temperature profile of the 
fiow [7 . The state of the fiow is sampled at 11 Hz by 
capturing raw intensity images [Fig. [l] (a) & (b)] using a 
12-bit digital camera; a background image of the fiuid be- 
low the onset of convection is then subtracted to remove 
optical nonuniformities (e.g., optical imperfections in the 
bottom plate). Digital Fourier filtering is then applied 
to remove high wavenumber components due to camera 
spatial noise. The data is then converted to binary values 
by thresholding at the median value of intensity. The re- 
sulting images [Fig. [l] (c) & (d)], where black represents 
hot upflow and white represents cold downfiow, are used 
as input to computation of homology. In what follows, 
X(n, e) represents the n^^ binary image in a time series 
at the reduced Rayleigh number e. Moreover, X^^^(n,e) 
and X^^^^(n, e) denote the hot fiow and cold fiow regions, 
respectively, of X(n, e). It is worth noting that the num- 
ber of pixels in X^^^(n, e) is equal to the number of pixels 
in X^^^^(n, e), i.e., hot and cold regions occupy the same 
area in the pattern. 

Homology theory provides a rigorous, systematic, and 
dimension independent method for characterizing geo- 
metric structures in X(n, e) using a few numbers. More 
specifically, the homology of a structure X in two dimen- 
sions (e.g., X^^*(n,e) or X^^^^(n,e)) is characterized by 
two nonnegative integers i — 0,1 called Betti num- 
bers, where /3o counts the number of connected compo- 
nents (pieces) of X, and Pi is equal the number of holes in 
X. A concrete illustration of this classification applied 
to convection patterns is given in Fig. [5] (For a more 
detailed discussion see jH |9|.) The computations were 
done using the package CHomP [9l [10] . Typically it takes 
about 10 seconds to reduce each shadowgraph image (of 
size 515 X 650 pixels) to binary form and to compute its 
homology. 

Parameter Distinction. — The sequence of Betti num- 
bers can be used to clearly distinguish different complex 
states of SDC (Fig.|3|. For the Rayleigh numbers e ^ 1.0 
and e ^ 2.0 we calculated sequences of Betti numbers 
/3f ^^(n, e) and /?f^^^(n, e), n = 0, . . . , 50, 000. Figure [s] 
shows that the time average of the Betti numbers are 
different at different e. 

The difference in the mean values of the Betti numbers 
(Fig.|3| refiects the instability mechanisms operating dur- 
ing the evolution of the complex spatiotemporal pattern 
in SDC. The dynamically significant events for the evo- 
lution occur in regions of local compression or dilatation 
of the roll structure. The compression leads to the merg- 
ing of neighboring rolls while the dilatation results in the 
formation of a new roll in the pattern. Such local events 
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FIG. 2: (Color online) Betti numbers for the region 
shown in Fig. [l](c). Image (a) shows the /3q^'^ = 34 distinct 
components, with neighboring components distinguished by 
different colors. Image (b) indicates the boundaries of the 
/3f = 13 holes. The only components, out of the /3§°^^ = 
38, of X^'^^^ (white region) that do not touch the boundary 
are the ones enclosed by the pi'''^ = 13 loops in X^''*. So 

the number of components of X^°^^ touching the boundary is 

^coid _ phot ^ 25. 



change the topology of the pattern. For instance, the 
merging of two rolls reduces the number of distinct rolls 
reducing Pq by one. A pattern increasingly dominated 
by such instabilities would hence show an overall reduc- 
tion in Po as different components link together locally 
at distinct locations. As the number of distinct rolls re- 
duces, such linkages are often self-intersections of a roll. 
This manifests as an increase in the number of loops {Pi) 
in the pattern. The mechanisms characterizing the sec- 
ondary instabilities for ideal straight rolls [TJ have a sim- 
ilar structure and are known as the skew-varicose (com- 
pression) and the cross-roll (dilatation) instabilities. In 
spiral defect chaos such events are seen to be operational 
in regions localized by the curvature of the rolls. 

Entropy. — At different parameter values, the state of 
SDC converges to different attr actors. The signature of 
convergence to the attractor is the evolution of the sys- 
tem in a bounded neighborhood in the space spanned by 
the Betti numbers. The set of Betti numbers describes 
the different patterns of rolls attained by the system. A 
complete description of the homological configuration of 
the state is determined by the set of Betti numbers de- 
scribing the hot and the cold regions. 

We calculate the probability of being in a particular 
configuration from a long time series of Betti numbers 
for the evolution of the state while on the attractor. The 
probability distribution enables the computation of an 
entropy to distinguish between different SDC attractors. 
The entropy is defined as 

i 

where the index i span the different states attained by 
the system as characterized by the Betti numbers. We 
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FIG. 3: Time series of Betti numbers showing every 400th 
data point of the time evolution of the number of distinct 
components and holes in the regions X^"^ (closed, black sym- 
bols) and X^"^^ (open, white symbols). The complete time 
series plots are qualitatively similar to the ones shown, the 
data is partially sampled only for better visualization. Circles 
correspond to f3o and diamonds to f3i . The horizontal lines in- 
dicate the mean values ^i{e) for the corresponding time series. 
Image (a) is for e ^ 1.0 and (b) for e ^ 2.0. 



nents do not interact, the entropy reduces as a result of 
no fluctuations in the Betti numbers associated with de- 
fect formation. The entropy increases for attractors that 
show evolving topologies mediated through defects. We 
flnd that the entropy of the system increases from 8.927 
to 9.631 as e is increased from e ~ 1.0 to e ~ 2.0. 

Stochastic Evolution. — The evolution of complex ge- 
ometries may be distinguished as being chaotic or 
stochastic from the time series of Betti numbers. The se- 
quence of Betti numbers has been used to uncover global 
chaotic evolution through the computation of the largest 
Lyapunov exponent in numerical simulations of reaction- 
diffusion systems previously [8^. In our experiments on 
SDC we have been unsuccessful in extracting Lyapunov 
exponents using similar techniques. In the case of SDC 
a mechanism describing locally chaotic islands driving 
the complex dynamics has been proposed [5^ It is 
likely that the spatiotemporally localized nature of such 
instabilities in SDC decorrelates at a scale that is not 
effectively captured by the time series of Betti numbers; 
a succession of these local events interspersed across the 
system may cause effectively stochastic evolution for the 
global geometric structure attained at the attractor. To 
flrst order we find it likely that the dynamics of fluctu- 
ations in the Betti numbers may be primarily stochastic 
in nature, as it is also suggested by the auto-correlation 
of the time series (Fig. |4|. 



2> 0.5 

i- 

o 

o 

o 

I 





8-1.0 





J 


£-2.0 





-5 5 



denote by pi the probability of being in a particular state 
described by /^f^S/^g^^'^ and /^f^^'^. The probability 
is obtained by trivially counting the number of distinct 
configurations attained over a large (ergodic) period of 
time and normalizing by the total number of possible 
configurations. While the entropy is computed over the 
states attained by the system in a finite amount of time, 
its value is seen to asymptote to a constant for long time 
series as the probability distribution of the Betti numbers 
at the attractor converges. In states where the compo- 



FIG. 4: The autocorrelation for the time series of the num- 
ber of holes (pi) in a state of SDC. The spiked nature of the 
function suggests that the fluctuations are primarily stochas- 
tic with little correlation. The autocorrelation for /3o has a 
similar spiked shape. 

Symmetry Breaking. — The homology of the states ex- 
hibited by SDC uncovers a break in the symmetry be- 
tween upfiows and downfiows. For an ideal Rayleigh- 
Benard problem, as described by the Boussinesq equa- 
tions [T2|, the inversion of the sign of the velocity (and 
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time) is also a solution to the evolution equations. This 
symmetry of the flow (the Boussinesq symmetry) in an 
ideal Rayleigh-Benard problem suggests that the statis- 
tical properties of the patterns would be the same for 
X^^^ and X^^^^] in particular we should obtain /3f^*(e) = 
^co^d(^) for z = 0,1. As Fig.H indicates, in our experi- 
ments we found that that the mean values of the Betti 
numbers clearly distinguish between fluid regions com- 
prised of hot and cold flows. Furthermore, this asymme- 
try is enhanced with an increase in the Rayleigh number. 
We suspect that the pattern homology serves as a sen- 
sitive detector of non-Boussinesq effects that are present 
due to the variation in physical properties of the fluid be- 
tween the hot bottom layer and the cool top layer. The 
strength of non-Boussinesq effects in experiments can be 
estimated by a dimensionless parameter Q [11, 13 com- 
puted perturbatively at the primary instability for con- 
vection. In our experiments with CO2, we find this pa- 
rameter to be equal 0.54 near the onset of convection 
— thus indicating strong, 0(1), non-Boussinesq effects. 
A similar computation at e ^ 1.0 and e ~ 2.0 yield Q 
equaling 1.04 and 1.49 respectively. 

Boundary Effects. — The pattern homology can also 
provide a well-defined way to separate boundary-driven 
effects from bulk phenomena in pattern forming systems. 
Typically, separating boundary behaviors from bulk ef- 
fects is done by setting a cutoff based on pattern correla- 
tion lengths. Using homology, we distinguish regions of a 
pattern as being part of the bulk if they are isolated from 
the boundary. These bulk regions of a given component 
(say, for example, isolated hot regions) are easily counted 
by recognizing they comprise the interior of holes of the 
other component (in this example, cold holes, as char- 
acterized by P^^^^). Thus, the number of components of 
Xhot connected to the boundary is p^^^^ ^ := - /3f 
(see Fig. H). Similarly for X^^^^, p^fy^^^ := - p^^\ 

One notices an asymmetry between the two components 
in this measure (Fig. |5| with a smaller number of com- 
ponents connected to the boundary as e increases for hot 
rolls as opposed to being almost the same for cold rolls. 
This also reflects in the number of components in the 
bulk, Po — Pbdy,o^ where the e dependence is primarily 
seen in the cold rolls. 

In summary, we have presented a robust topological 
characterization of experimental data from a spatiotem- 
poral complex dynamical system. This technique gives 
one the ability to map the gross global configuration of a 
state to a few dimensionless numbers. We found stochas- 
tic (and ergodic) global topological dynamics for the sys- 
tem of spiral defect chaos where the underlying micro- 
scopic dynamics is deterministic and previously observed 
to be locally chaotic. We also found asymmetric topologi- 
cal configurations between hot and cold regions for which 
we are unable to provide a physical mechanism. A more 
complete description of the system would require cou- 
pling such scale-independent measures with other mea- 




12 3 4 



" xio' 

FIG. 5: The number of components connected to the bound- 
ary for e ^ 1.0 (a) and e ^ 2.0 (b). As in Fig. ISlclosed circles 
represent and open circles represent X^° ^ . A sample 

from every 400th data point from the complete time series is 
displayed in the above figures. 



sures that depend on the length-scales of the patterns, 
e.g. the mean length of a roll between defects. Addition- 
ally, it would be interesting to compare the correlation 
of local characterizations such as Lyapunov exponents or 
wave numbers with variations in quantitative topologi- 
cal measures. Understanding the interplay between local 
and global dynamics is a core issue in the study of com- 
plex systems. 
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